IAES International Journal of Robotics and Automation (IJRA) 
Vol. 11, No. 1, March 2022, pp. 43~61 
ISSN: 2722-2586, DOI: 10.11591/ijra.v11i1.pp43-61 og 43 


An effective approach to enhance the balancing control in 
bycycorobot using the soft computing techniques 


Aswant Kumar Sharma, Dhanesh Kumar Sambariya 
Department of Electrical Engineering, Rajasthan Technical University, Rajasthan, India 


Article Info ABSTRACT 

Article history: The balancing and control of bycycorobot is a challenging task. The pre- 
. specified controller available in the literature for balancing has been reduced 

Received Nov 8, 2021 with novel optimization to improve the effectiveness of balancing, 

Revised Dec 25, 2021 uncertainty, and the complexity of the complete system. The novel Harris 

Accepted Dec 30, 2021 hawk optimization (HHO) which is based on the hunting behavior of the 


hawk has been utilized to improve the balancing of the bycycorobot. The 

paper proposes the decreased order controller of a pre-specified controller 
Keywords: for a bycycorobot. The obtained controller response with bycycorobot in the 
complete closed loop is analyzed, and the best performance is compared 
with the reduced order controller available in the literature. The comparison 
is based on the response indices and response characteristics. 


Balancing control 
Bycycorobot 
Harris hawk optimization 


Reduced controller 
Reduced-order modelling This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Aswant Kumar Sharma 

Department of Electrical Engineering, Rajasthan Technical University 
Akelgarh, Kota, Rajasthan 324010, India 

Email: aswantksharma @ gmail.com 


1. INTRODUCTION 

As the first kind of personal mobility vehicle ever invented, bicycles play an important role in the 
history of transport. Bicycles are lightweight and are solely propelled by human power. They can contribute 
to the reduction of traffic congestion and air pollution in urban areas. Bicycles are classified as single-track 
vehicles and display interesting dynamic behavior. The intriguing behavior of bicycles poses challenging 
problems in modeling and control, which have attracted attention from the automatic control research 
community [1]. The balancing control problem of the bycycorobots is not new in the field of robotics. The 
bycycorobot is in great demand due to its basic construction, simpler dynamics, and applicability in wide 
sectors such as transportation, security, search and rescue, and labor reduction. These robots are inherently 
unstable and susceptible to external disturbances. Robust controls approached are required for appropriate 
and smooth balancing control and movement of such robots [2]. The stability of bicycles is an issue that has 
scarcely been understood in dynamics and is a commonly ignored problem. After two centuries of debate and 
unfinished modeling, new researches are paving the groundwork for increasing studies regarding this issue 
[3]. Thereby, some researchers have appeared thanks to new dynamics models and the computational power 
available nowadays. Their purpose is to propose and test different alternatives to stabilize a bicycle. Despite 
their common goal of making a bycycrobot stable, they propose a wide variety of applications for their 
stabilized bycycrobot [4]. 

The bicycle can be balanced using the rotation of the front wheel in the direction of lean, which 
changes the tire contact point with the ground in the same direction and is similar to balancing an inverted 
pendulum. In addition, the centrifugal force contributes to balance because of the circular motion. Moreover, 
uncontrolled bicycles can balance themselves within certain velocity ranges, which depend on different 
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bicycle-design parameters [5]. For example, gyroscope, perturbation stationary bycycrobot system, inverted 
pendulum nonlinear controllers, robust two-stage active disturbance rejection control, fuzzy sliding-mode 
underactuated control, proportional derivative (PD) controllers, and proportional integral derivative (PID) 
controllers have been presented in distinctive researches [6]-[11]. The PID controller has three parameters to 
be tuned suitable for bycycrobot because of their degree of freedom-one for tilt, other for speed, and subject 
upon the applicability of the yaw. Numerous techniques are given in the literature for balancing control of 
bycycorobot like a flywheel, mass, and steering balancing [12]-[14]. Although much effort has been done to 
understand this fascinating phenomenon, the studies now have been restricted to a single situation in which 
the bicycle rides on a horizontal surface [15]. The adaptive integral terminal sliding mode was used to 
provide robust control of the reaction wheel of bycycorobot [16]. The proposed adaptive integral terminal 
sliding mode (AITSM) control scheme can not only stabilize the bicycle robot and reject external 
disturbances generated by uncertainties and un-modelled dynamics, but also eliminate the need of the 
required bound information in the control law via the designed adaptive laws. The experimental results verify 
the excellent performance of the proposed control scheme in terms of strong robustness, fast error 
convergence in comparison with other control schemes. The “scissor-pair control moment gyroscope 
technique” [17] for active balancing of the unmanned bicycle system. A miniature bicycle system was 
constructed, and its 3D solid model was used to acquire the dynamics model’s parameter values. There are 
two types of disturbances that might cause a bicycle with an unstable equilibrium point to become unstable 
impulsive external disturbances and static continues disturbances. The experimental result depicts the 
performance of active balancing for a bicycle equipped with a scissor-pair control moment gyroscope (CMG) 
in the presence of such disturbances. 

A bicycle robot balancing with steering and a non-inverted pendulum [18]. Euler-Lagrange equation 
of motion was used to create a nonlinear mechanical model of the bycycorobot moving at a constant forward 
speed, which was then linearized around the upright position. The pendulum and steering motor models were 
also included in the linearized mechanical model to obtain the complicated voltage-input system represented 
by the state-space model. The electro-mechanical parts of the robot were designed considering the controller 
for stabilizing the system. The simulation output revealed that the system with the pendulum-emphasized 
gain matrix outperformed the system with the steering-emphasized gain matrix in terms of learning and path- 
tracking ranges. The simulation results were, likewise, in line with the real outcome on the straight road and 
turning points. 

One of the major problems in the development of bycycorobot is balancing at modest forward 
speeds. Therefore, Tofigh et al. [19] demonstrated a novel gyrostabilizer with a dual-flywheel configuration 
that can give any required gyroscopic roll moment. Extensive simulation comparison has been being carried 
out to assess the approach and demonstrate the impact of incorporating the new actuator. The result in [19] 
showed that the proposal provides higher performance while keeping the control effort within the capabilities 
of standard actuators. A structure-specified mixed “H2/Hoo” controller was designed using particle swarm 
optimization (PSO) for control balancing of bycycorobot. The unstable system with many sources of 
uncertainty due to un-model dynamics, parameter variations, and external disturbances [20]. The design 
process, on the other hand, frequently generates a complicated and non-convex optimization issue that is 
difficult to address using traditional optimization approaches. PSO was a recently used meta-heuristic search 
approach for solving non-convex and multi-objective optimization problems. PSO was utilized in the 
approach to find parameters of a structure-specified controller that satisfy the combined H2/Ho performance 
index. The simulation and experimental findings demonstrate the proposed controller’s resilience when 
compared to the traditional PD controller, as well as the suggested algorithms’ efficiency when compared to 
the genetic algorithm (GA). Sikandar and Prasad [21] suggested combining reduced order modeling and a 
cuckoo search algorithm to develop a pre-specified structural controller for balancing control of two-wheeled 
mobile robot. A lower order pre-specified structural controller was created in this study, and it was found to 
be efficient enough to deal with unknown dynamics. The undetermined parameters of the proposed reduced 
order controllers are determined using the cuckoo search method. The suggested controllers’ results s was 
compared to those of GA, PSO, Schur analysis, balance truncation, model truncation and traditional PD 
controller. The proposed controller was found to have superior performance when compared to other 
controllers. In terms of time response requirements and performance indices such as integral of square error 
(ISE), time-weighted integral of absolute error (ITAE), integral of absolute error (IAE), the performance of 
higher and lower controllers was also investigated using a perturbed two-wheeled mobile robot. 

Controlling the equilibrium of a bycycorobot is a particularly tough task. Many balance control 
methods for two-wheeled bicycles have been suggested including non-linear control, compensator design 
utilizing the roots’ locus techniques, PD control, PID control, and fuzzy control. However, these control 
algorithms are not robust; the two-wheeled bicycle model cannot bear varying weight loads and cannot bear 
varying weight loads and cannot operate in noisy situations. As a result, a robust control method is 
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appropriate for the real-world applications of the two-wheeled bicycle model for this, composite method 
H2/Hoo control is presented, which is a sophisticated approach for constructing durable and optimum 
controllers for systems with unknown sources. However, it is widely understood that the structure-specified 
mixed H2/Hoo controller design typically creates a complicated and non-convex optimization issue that is 
difficult to address using traditional optimization methods. The optimization approaches using GA, PSO, and 
Cuckoo search have been suggested to tackle this problem. These techniques have offered to balance but the 
improvement of balancing also has the possibility in bycycorobot against uncertainty and external disruption. 
The reason of unmolded dynamics, parameter changes, and external disturbances complicates the system and 
necessitates a strong controller. As a result, a novel approach is required to enhance the balance control in 
bycycorobot. 

The model order reduction methodology has been applied in order to obtain the reduced order 
controller. The MM techniques [22], continued fraction [23], mixed methods [24]-[28] are traditional 
methods tested and failed. The main objective of this manuscript is to obtain the effectively reduced order 
controller with bycycorobot response and enhance the balance control of it under uncertainty and external 
disturbances. For this, Harris hawk optimization (HHO) is having been selected due to aggressive and swarm 
behavior to design the controller efficiently. The approach represents the systems’ model uncertainty as 
multiplicative uncertainty, and the system is considered to be influenced by external disturbances. HHO has 
been utilized for search parameters of and admissible structure-specified controller that minimizes the ISE 
while being subjected to robust stability constraints (Hoo norm) against model uncertainty and external 
disturbances. The suggested method is used to regulate the balance of a bicycle robot equipped with a 
gyroscopic stabilizer, known as bycycorobot. Using the Lagrange technique, a simplified dynamics model of 
a bycycorobot is constructed by disregarding forces caused by forwarding movement and steering. To 
evaluate the performance proposed system, it was compared with other techniques available in the literature. 
The proposed system is implemented using MATLAB software. The proposed algorithm-based controller 
will be effectively enhanced the performance of the balancing control in bycycorobot which is an unstable 
system with un-modeled dynamics, parameter changes, and external disturbances as sources of uncertainty. 


2. HARRIS HAWK OPTIMIZATION (HHO) 

HHO is based on the studies of hawk behavior usually in the period of hunting. The study has been 
carried out by Louis Lefebvre. The mathematical implementation in the engineering designing field using the 
HHO is carried out [29]. The behavior of hunting and chasing patterns for the capture of prey in nature is 
known as a surprise pounce. The searching of prey is a task done by the predator using the highest point of 
the area such as standing on top of trees or flying in the sky. The attack of the hawk on prey is called a 
pounce. As the prey is spotted another member is informed by visual displaying or vocalization. The HHO is 
divided into three-phase naming exploration, the transition from exploration to exploitation, and the 
exploitation phase. The exploitation stage is separated into four stages namely soft besiege, hard besiege, soft 
besiege with advanced quick dives, and hard besiege with progressive speedy dives. Figure | is showing the 
flow chart of HHO. 


2.1. The exploration phases 

To start this phase, the Hawk reaches the peak of tree/pole/top of the hill to trace the prey and also 
consider the other of Hawks positions. The situation is q > 0.5 of branch on random giant trees for the 
situation of q < 0.5. The condition ids are modeled as (1). 


Xrand (t) = rılXrana (t) ~ 27r,X (t)|q 20.5 


Xprey(t) — Xm(t) — 74(LB + r,(UB,p))) q < 0.5 (1) 


x(t+1)= k 


Where X(t + 1) is the position vector of the hawk in the succeeding iteration t.Xprey(t) is the present 
position vector of hawks r4, T2, T3, r4 and q are the random number confidential (0,1) upgraded with iteration. 
LB is the lower bounds, and UB is the upper bounds of numbers X;gnq(t) represents the arbitrarily hawk 
represents the arbitrary hawk from the present population Xm is the average position of the current population 
is the average position of the current population of hawks. The primary rule creates solutions based on a 
random position. In the second rule of (1), the variance between the best positions and the average location of 
the group plus an arbitrary climbed factor depends on the number of variables. The scaling factor r} increases 
the random nature of regulation once 7 adjacent value to 1 adjacent value to 1 and comparable distribution 
designs. Random factor scaling coefficients increase pattern diversification and explore various feature 
regions. The rules for buildings are capable of mimicking the actions of a hawk. The hawk’s average location 
is obtained using (2). 
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Where X,,,(t) is obtaining by (2) X;(t) designates the position of an individual hawk in iteration f and N 
signifies the number of hawks. 


2.2. Conversion from exploration to exploitation 

The exploration to exploitation changes between exploitation performances founded on the 
absconding energy of the prey. The energy of a prey reduces throughout the escaping. The energy of the prey 
is modeled as (3). 


E = 2E)(1--) (3) 


Where E designates the absconding energy of prey, T is the maximum number of iterations, and Eo is the 
initial state of energy. 


2.3. Exploitation phase 

The process begins with the surprise, and the imagined prey of the previous stage is hostile. The 
preys try to escape. The probability of fleeing from the prey is (r < 0.5) or not to escape efficaciously 
(r = 0.5). The hawk executes rough or soft besieges concerning prey activity to capture the prey. Based on 
the vitality of the prey, the hawk encircles around the beast in various ways. The hawk gets closer to the 
desired prey to maximize its odds of cooperating in killing the rabbit. The gentle assault begins, and the 
rough assault takes place. 


2.4. Besiege occurs 
Besiege is the process at the time of capturing prey. It is divided into soft besiege and hard besiege 
with progressive dives of each respectively. 


2.5. Soft besiege 

The prey has energy and tries to escape using random confusing jumps. The value for escaping 
energy must be r > 0.5 and E = 0.5. If the values are below as stated, the prey is unable to jump. Hawk 
encircles prey gently to make it more tired and achieve the surprise dive. This conduct is modeled by 
subsequent rules represented in (4) and (5). 


X(t + 1) = AX(t) — E|J * Xprey(t) — X(t)| (4) 

AX(t) = Xprey (t) — X(t) (5) 
2.6. Hard besiege 

The prey is exhausted and has less energy when values r > 0.5 and E > 0.5. The hawk barely 
encloses the intended prey and finally achieves the shock pounce. The present locations are updated as (6). 

X(t + 1) = Xprey(t) B E|AX(t)| (6) 

Dive is founded on the LF-based designs using the law represented in (7). 

Z=Y+SxXLF(D) (7) 


Where D represents the dimension problem and S represents a random vector by size 1 x D and LF is the 
levy fight function, and calculated as (8). 


x r(atp)xsin(™4) i 
LF (x) = 0.01 x2,0= = (8) 
|v|2 r(22)xpxall2 


Where, u,v are random values inside (0,1), 8 is a constant set to 1.5. The last tactic for apprising the 
locations of hawks. The soft besiege stage can be achieved and given in (9). 


IAES Int J Rob & Autom, Vol. 11, No. 1, March 2022: 43-61 


IAES Int J Rob & Autom ISSN: 2722-2586 


Y if F (Y) < F(xX() 


IEEE a if F(Z) < F(X(t)) 


The Y and Z are obtained using the (8) and (9). 


Define HHO parameter: 
Population size; max. iteration 
optimization dimension 

¥ 


Construct a population 
(randomly) 


Set itr=1 


—— 


calculate objective function(ISE) 
choose best solution 


Update: 
E = 2rand(0,1)-1,J =1(1 -rand (0,1)) 
and E using defined Egns. 


Update the solution 
by Eq. (1) 


Update the solution 
by Eq. (2) 


Update the 
A solution by Eq. (3) 


Update 
t=t-+1 


t Update the solution 
by Eq. (.5) 


Update the 


— 


Report best solution and 
its objective function 


Figure 1. Flow chart of HHO 


2.7. Hard besiege with progressive dives 


(9) 


The |E| < 0.5 and r < 0.5. To escape and hard besiege is built earlier the surprise pounce to catch 
and kill the prey. The condition on the prey side is comparable to that of soft besiege except this time. The 


hawk seeks to reduce the difference between their regular position and the fleeing target. 


_ (Vif F(Y) < F(X) 
Ae Lye z if F(Z) < F(X(t)) 
The Y and Z are gained by (11) and (12). 
Y= Xprey(t) = E|J X Korey (t) = Xm(t)| 


Z =Y +S xLF(D) 


(10) 


(1) 


(12) 
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3. MATHEMATICAL MODEL OF BYCYCOROBOT SYSTEM 

The bycycorobot is shown in Figure 2. Figure 2(a) is showing bicycle representation of bycycorobot 
for side view, Figure 2(b) is showing bicycle representation of bycycorobot for front view, Figure 2(c) is 
showing bicycle representation of bycycorobot for flywheel side view, Figure 2(d) is showing bicycle 
representation of bycycorobot for flywheel front view. Wheels are mounted on a different axis. The main task 
the goal of the robot is to move without falling in any direction, with or without limited load. The major issue 
of the taken task is of balancing because of its unstable nature and various uncertainties. Different 
controllers/ algorithms are available to resolve the problem. Bycycorobot dynamics is derived using the 
Langrage (13). 


dt 


d ga OKE OPE 
əfi 


Cu oe = 


Where KE denotes the total kinetic energy, PE represents the total potential energy of the system, fiis 
generalized coordinate, and EF;denotes the external forces. The PE and KE relation are (14)-(17) 


Bycycorobot dynamics are derived using the Langrage (14). 


PE = m, gH, cos a + mpghp cos a (14) 
1 x $ 4 1 ` 1 5 PARS, PE 
KE =>M,(&? H?) + 5M, (&? he) + zl, P+ 5 mB? + Ip (á sin B)? + Im(a sin B)*| (15) 


fi=a (16) 


_ [m,H? + mph; + 1, + 

I, sin? B + Im cos? B 
+2 sin B cos B (I; — Im )aB — 
g(m,H, + mghr) sina 


=IpwB cos Bp (17) 


Where, M, and My represents the weight of bycycorobot and flywheel respectively. a represents the lean 
angle, & denotes angular velocity in Z-axis, £ represents the angle of flywheel along Z1 axis, Å is the angular 


velocity of the flywheel along X; axis. H, is the height of center of gravity of bycycorobot, hy represents the 


height of flywheel center of gravity. Im, Ir, and I, are flywheel radial movement of inertial, flywheel polar 


movement of inertial and robot movement of inertia respectively. The value of f; = f and the equation 
obtained in (18). 


Blin — @? (Ip — Im) sin B cos B = Tm — Ipwa cos B — Vin B (18) 


The chain transmission of bycycorobot and DC Motor dynamics is assumed to be 5:1 and the 
relation obtained is given in (19) and (20). 


Tm = 5Kmi (19) 
U=LĒ+Ri+ KÊ (20) 
Where in (19), Km represents the torque constant of motor and in (20) Ke is back e, m, f. constant, R and L 


are armature resistance and inductance of motor. The equation (19) is substituted in (18) and linearizing (17) 
and (18) around the equilibrium point, the relation achieved as (21) and (22). 


$ Mh, \. , 
&[M,H? + Mph; + 1, + Im] -g (sick) å- lw = 0 (21) 
Blin + Ipod + VinB — 5Kmi = 0 (22) 
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mi; għ, cos a 


m gh, osa 


(b) 


(c) (d) 


Figure 2. Bicycle representation of bycycorobot for (a) side view, (b) front view, (c) flywheel side view, and 
(d) flywheel front view 


The consideration of the equation x = [a & B i] and y = q along with u = U. The state space 
equation of the system form is given in (23). 


x = Ax + Bu 
k = Cx + Du (23) 
The values of the fourth order state space input part, output part is represented as (24)-(27). 
0 1 0 0 
g(mrhr+mphf) Ifo 
EE A 0 EEE TA 
mrhr+mfhý+ir+im mrhy+mfhý+lr+Im 
A= Ifo Vm 5Km (24) 
0 STA —m ohm 
Im Im Im 
0 0 ae -2 
L L 
a 
B= [o 0 0 ] (25) 
C=[1 0 0 0] (26) 
D = [0] (27) 


In most cases, a flywheel is used to balance the generated torque with respect to gravity. The 
parameter consideration of bycycorobot is listed in Table 1. 
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Table 1. Parameter values 


Parameters Values Parameters Values 
Mass of the robot (m,.) 8.1 kg Height of robot centre of gravity (H,) 0.86 m 
Height of robot centre of gravity (H) 0.86 m Robot moment of inertia 27.548 kgm? 
Robot moment of inertia 27.548 kgm? Mass of flywheel (m,) 43.1 kg 
Flywheel polar moment of inertia 0.215926 kgm? Height of flywheel centre of gravity 0.8m 
Torque constant of the motor (Km) 0.119 Nm/A Flywheel radial moment of inertia 0.112304 kgm? 
Flywheel speed (w) 157.08 rad/s DC motor viscosity coefficient (V,,) 0.000253 kgm?/s 
Inductance of motor (L) 0.0006 H Back e.m.f constant of motor (K,) 0.1184 Vs 
Mass of the robot (m,.) 8.1 kg Armature resistance of motor (R) 0.41 


The state space into transfer function. The obtained transfer function of the two wheeled mobile 
robot is given as (28). 


G(s) = ADAE EE ee wk (28) 


U(s)  s44+683.353+120852+109700s—6949 


The output lean angle and the input voltage to the DC motor that controls the flywheel control axis. 
Assuming the following two cases: i) Case-1: the additional 10 Kg load and decrease the speed of the 
flywheel up-to 147 rad/sec. The transfer function from this condition is as (29); ii) Case-2: in this case, more 
10 Kg load is added speed of flywheel is increased to 167 rad/sec. the transfer function is as (30). 


3784 
s*4+683.3534+1162s2+78290s—6857 


G(s) = (29) 


4299 
S*4+683.353+1197s2+102300s—6857 


G(s) = (30) 


Bycycorobot including both the special case of the system represents instability in Figure 3. The 
balancing control of the system is given by Thanh and Parnichkun in [20] with a controller based on particle 
swarm optimization using specifically mixed H2/Ho controller. The transfer function of the controller is as 
(31). 


1275s°+8.695X10°s*+5.151X10°s3 
G (s) = +1.359x1085?+2.435x1075+1.091x10° 
(Hoo)e S6+715.78°+2.355X10454+42.789X10953 
+3.802X10°s2+6.519X10°s+2.872X10* 


(31) 


4. DESIGN OF FIRST ORDER CONTROLLER USING HHO TECHNIQUE 
The structure of first order controller is obtained using the HHO by minimizing the ISE. The 
unknown reduced order model of first order is represented as (32). 


Ny 


Gc-Pr opsoed (s) = (32) 


The unknown value of the proposed first order controller is optimized using the Harris hawk 
algorithm. The important parameters lower bound (LB) is [150 1], upper bound (UB) is [200 160] number of 
iterations 100, search agents or hawks is 10. The obtained first order controller after 100 iterations in 
353.148409 s is as (33). 


191.1806 


1+5.0382 (3) 


Gc-Pr opsoed (s) = 


The numerator, denominator parameter of first order controller optimized using the Harris hawk 
optimization and obtained result shown with iteration graph and response of proposed first order controller 
with sixth order and first order controller from literature given in Figure 4. 1* order controller with 6" order 
H-œ and 1“ order controller from literature for: i) Numerator as shown in Figure 4(a); ii) Denominator 
parameter as shown in Figure 4(b); iii) Fitness function as shown in Figure 4(c); and iv) Step response as 
shown in Figure 4(d). 
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Figure 3. Step response of the bycycorobot (original, case-1 and case-2) without controller 


« D} = 5.0382 
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(c) (d) 


Figure 4. Proposed 1“ order controller with 6" order H-o and 1“ order controller from literature for 
(a) numerator, (b) denominator parameter, (c) fitness function, and (d) step response 


The H-co 6" order controller response, shown in the blue step response and it contains some 
uncertainties and disturbance in the period of 0-1.5 s. The proposed first order controller using the HHO 
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follows H- controller closely and the uncertainties and disturbance are also removed. This indicates that in 
closed loop with bycycorobot in all three cases, it may perform with a better result. For Justification the 
response is also compared with the first order controller using GA and PSO [20] given as (34) and (35). 


197.33 

Gc-Ga = S+4.91 (34) 
196 

Gc-pso = 344.58 (35) 


Recently the first order controller with improved performance using the cuckoo search optimization 
in [21]. 


150 


Gc-csa =a 367 (36) 
The first order controller with improved performance using PSO, the structure specified Hoo loop 
shaping controller in [30] and is given as (37). 


5 
Gopo == (37) 


S+4.63 


The proposed controller needs to be tested and the response of it is compared with the controller 
available in the literature. 


4.1. First order controller closed loop analysis with bycycorobot 

The proposed first order controller in closed loop with the bycycorobot including the special case-1 
and case-2 has been analyzed in this part. The performance has been analyzed on the basis of step response 
characteristics and response indices error in Figure 4. The step response of the 1* order controller using HHO 
with bycycorobot in a complete closed loop is represented in Figures 5-7 respectively. 

The performance in terms of step response characteristics and response indices error of proposed 
first order controller with the system in closed loop with all cases available in Table 2. The proposed 
controller response is better than the first order controller using GA and PSO [20], Cuckoo search 
optimization in [21], the structure specified Hoo loop shaping controller using PSO [30]. This can also be 
verified from the step response plot in Figure 5, Figure 6 for case-1, and Figure 7 for case-2. 


Table 2. Closed-loop performance of proposed 1“ order controllers with bycycorobot and its case-1, case-2 
Response indices error 


System with controllers 


ISE IAE ITAE ITSE 
OEA 0.008372 0.1553 0.4408 0.007561 
Case 1 0.008372 0.1553 0.4408 0.007561 
Case 2 0.007634 0.1556 0.4703 0.007127 
Gc-ça [20] 0.01018 0.225 1.461 0.01385 
Case 1 0.0118 0.2547 1.766 0.01629 
Case 2 0.009557 0.2346 1.624 0.01453 
Gc-pso [20] 0.02176 0.8046 31.27 0.2131 
Case 1 0.01196 0.255 1.753 0.01634 
Case 2 0.009659 0.2349 1.613 0.01455 
Gc-pso1 [30] 0.03775 1.343 61.73 0.7739 
Case 1 0.03299 0.5856 7.03 0.1223 
Case 2 0.02768 0.5173 6.036 0.0926 
Gc-csa [21] 0.0136 0.1547 0.2166 0.00881 
Case 1 0.01759 0.2749 2.342 0.02219 
Case 2 0.01312 0.2366 1.964 0.01599 


The data available in Table 2 concludes that the proposed controller is efficient and effective in all 
two cases and performs better. The controller analysis extended to the second order and its design using the 
HHO. 
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Figure 5. Step response characteristics of the Figure 6. Step response characteristics of the 
proposed 1* order controller in closed loop and first proposed first order controller in closed loop case-1 
order controller from literature [20]—[22] and first order controller from literature [2]—[4] 
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Figure 7. Step response characteristics of the proposed first order controller in closed loop case-2 and first 
order controller from the literature [2]—[4] 


5. DESIGN OF SECOND ORDER CONTROLLER USING THE HARRIS HAWK 
OPTIMIZATION 
The design of unknown second order reduced controller based on the structure of bycycorobot is 


represented as (38). 


NysS+N2 (38) 


Gc- s) = 
Cc Pr oposed2( ) s2+D1S+D2 


The response analysis of the system with higher order controller and reduced controller from the 
literature the lower bound selected as Ib = [230 950 8 20]; and upper bound selected as ub = [300 1000 10 
30] with N,, Nz, D, and D, are unknown four parameters dimensions which are to be optimized using an 
algorithm so dimension is 4. The obtained parameter using the Harris hawk algorithm is represented in 
Figure 8, where Figure 8(a) is showing optimized numerator (Nı), Figure 8(b) is showing optimized 
numerator (N2), Figure 8(c) is showing denominator parameter values (Dı), and Figure 8(d) is showing 


An effective approach to enhance the balancing control in bycycorobot using ... (Aswant Kumar Sharma) 


54 o ISSN: 2722-2586 


denominator parameter values (D2). The obtained second order using the HHO algorithm is represented as 
(39). 


G OE 292.22475+999.1637 i 
C—Pr oposed2 ~~ 5249.83635+26.3188 


Figure 9 is showing the iteration graph of the controller with the fitness value and step response the 
higher order Hoo controller with the proposed reduced order, reduced order available in the literature is given 
in Figure 9(a). The step response in Error! Reference source not found.9(b), it is clear that the proposed 2™ 
order reduced controller response approximately same as that of higher order controller and improved 
response from the controller from cuckoo search as (40) [21], and the second order controller using the PSO 
as (41) [30]. The effectiveness of the proposed order controller is also analyzed with the system in closed 
loop. 


266.55+891.6 
a (40) 
S“4+9.1395+25.03 


Gc-cs(s) = 


129.75+499.6 
Gc_pso(s) = $2+6.3835S+16.18 (41) 
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Figure 8. Harris hawk algorithm with a number of iterations for (a) optimized numerator (N1), (b) optimized 
numerator (N2), (c) denominator parameter values (D1), and (d) denominator parameter values (D2) 
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Figure 9. Iteration graph for (a) fitness value and (b) step response characteristics of higher order with 
proposed 2" order controller and 2™ order controller from literature 


5.1. Second order controller closed loop analysis with bycycorobot 

The second order controller in closed loop with the bycycorobot stabilizes the system by removing 
uncertainties. The response of the controller with response characteristics and response indices error is given 
in Table 3. 

The response of the proposed second order approximately gives the same response as that of the 6 
order controller with bycycorobot. Figure 10 is showing step response characteristics of the proposed second 
order controller in closed loop and second order controller from literature, Figure 11 is showing step response 
characteristics of the proposed second order controller in closed loop and second order controller from 
literature case-1, and Figure 12 is showing step response characteristics of the proposed second order 
controller in closed loop and second order controller from literature case-2. 

The response is also justified with the result from the cuckoo search [21] and particle swarm 
optimization [30] It is clear from Table 3 that the proposed controller provides the comparable performance 
and the response indices IAE, ITAE and ITSE is better than the result from cuckoo search and particle swarm 
optimization. 
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Table3. Closed-loop performance of proposed second order controllers with bycycrobot and its case-1 and 


case-2 

System with controllers Response indices error 

ISE IAE ITAE ITSE 
Gc-Pr opsoed 0.004346 0.1222 0.4305 0.004686 
Case 1 0.003511 0.1185 0.4399 0.003884 
Case 2 0.003414 0.1116 0.4181 0.003781 
Gc-cs 21] 0.003304 0.1284 1.025 0.004934 
Case 1 0.00256 0.1303 1.349 0.005456 
Case 2 0.002491 0.123 1.162 0.00482 
Gc-pso 30] 0.01567 0.3326 3.918 0.0401 
Case 1 0.02063 0.4242 5.222 0.06798 
Case 2 0.0157 0.3607 4.471 0.05039 
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= = = System with 2" order by CSO 2016 


System with 2"4 order by PSO 
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Figure 10. Step response characteristics of the proposed second order controller in closed loop and second 
order controller from literature 
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Figure 11. Step response characteristics of the proposed second order controller in closed loop and second 


order controller from literature case-1 
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Figure 12. Step response characteristics of the proposed second order controller in closed loop and second 
order controller from literature case-2 


6. DESIGN OF THIRD ORDER CONTROLLER USING THE HARRIS HAWK OPTIMIZATION 
The third order controller with unknown is represented as (42). 


N1S?+N2S+N3 
s3+D1S?+D2S+D3 


Gc-pr oposed2 (s) = (42) 


The value of unknown controller is optimized by HHO using Ib = [1240 230 193400 30 388 4470] 
and ub = [1270 240 193900 35 400 7970] with iteration of 100 and number of search agents are 10. The ISE 
is minimized in order to obtain the reduced order. 

The unknown parameter of reduced order model is optimized and the obtained parameter is with 
respect to number of iterations is represented in Figure 13 and third order obtained as (43), Figure 13(a) is 
showing optimized numerator (N1), Figure 13(b) is showing denominator parameter (D1), Figure 13(c) is 
showing optimized numerator (N2), Figure 13(d) is showing denominator parameter (D2), Figure 13(e) is 
showing optimized numerator (N3), and Figure 13(f) is showing denominator parameter (D3). 


1259.294172s7+233.19765+193766.8479 
s3431.1049s2+394.9840S+5276.281587 


Gc-pr oposed2 (s) = (43) 


The third order reduced system using the Schur analysis, balanced truncation and model truncation 
given in [31], is represented as (44), (45), and (46). 


1275s7+234.85+1.993e5 


G s)= 44 

c-sa(S) $3433.7852+395s+5506 (44) 
1275s7+233.8s+1.992e5 

Gc-pr(s) = =a" aa OLDS AEGa CaGe: (45) 
53433.78524+395s+5499 
1057s7+226.5s+1.6358e5 

Gc-mr(s) = ays Gast PROP WEPAEDy (46) 
53427.99s2+395.9s+4521 


The proposed reduced order is analyzed with closed loop with the bycycorobot and equivalent 
reduced order from the literature [31]. The third order controller by cuckoo search algorithm [21], is 
represented as (47). 


1241s7+234.75+1.993e5 
53432.47524+395s+5274 


Gc-cs(s) = (47) 
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Figure 13. The Harris hawk algorithm with number of iterations for (a) Optimized numerator (N1), 
(b) Denominator parameter (D1), (c) Optimized numerator (N2), (d) Denominator parameter (D2), 
(e) Optimized numerator (N3), and (f) Denominator parameter (D3) 


6.1. Third order controller closed loop analysis with bycycorobot 

The proposed third order response has been analyzed with the system in closed-loop. The response 
of the system with the proposed controller in special cases is mentioned in Table 4 and compares with the 
response with the controller available in the literature [21], [31]. The third order proposed controller gives the 
approximate same response as compared to the controller present in literature [21], [31] and represents with 
step response characteristics in Figure 14, Figure 15 for case-1, and Figure 16 for case-2. The response 
indices IAE, ITAE and ITSE are minimum indicating the good approximation and superiority of the 
proposed reduced order justified. The approach has been verified with the pole location of the proposed first 
order, second order and third order controller with the result from the literature [20], [21], [30], [31]. Table 5 
is showing poles of higher order and reduced order controller. 
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Table 4. Step response characteristics and response indices error of proposed reduced controller in closed 
loop with system 


1.2 


Amplitude 
= = 
n oO 


> 
z 


0.2 


System with controllers 


Response indices error 


ISE IAE ITAE ITSE 
Go-gřtpyopsded 8.121e-5 0.03723 0.4609  0.00005925 
Case 1 0.0001205 0.04514 0.6064 0.001008 
Case 2 8.782e-5 0.04059 0.5273 0.0007553 
Gc-sa [21], [31] 7.472e-5 0.04219 0.6628 0.001112 
Case 1 0.0001703 0.05152 0.8902 0.001961 
Case 2 9.283e5 0.04738 0.7647 0.001472 
Ge-er [21], [31] 7.205e-5 0.04131 0.6507 0.001073 
Case 1 0.0001032 0.05042 0.8733 0.00189 
Case 2 8.945e-5 0.04639 0.7505 0.00142 
Gc-mr [21], [31] 0.0002247 0.04865 0.6664 0.00117 
Case 1 0.0005654 0.07986 0.9448 0.002514 
Case 2 0.0002393 0.05389 0.7687 0.001534 
Gc-cs [21], [31] 0.0003682 0.05164 0.2913 0.0007287 
Case 1 0.0004263 0.05369 0.3004 0.000834 
Case 2 0.0003882 0.055 0.3187 00.0008373 
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Figure 14. Step response characteristics of the 
proposed third order controller in closed loop and 
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proposed third order controller in closed loop and 
third order controller from literature case-1 


Step Response 


System with 6" order controller 


Proposed 34 order controller by HHO 


3°4 order by SA 

wrd order by BT 
——3"4 order by MT 
rd order by CSA 


Time (seconds) 


Figure 16. Step response characteristics of the proposed third order controller in closed loop and third order 
controller from literature, case-2 
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Table 5. Poles of higher order and reduced order controller 


Controller Poles Controller Poles 
6" Order -681.74+0.0000i; -26.71+0.0000i; -3.54+13.92i; Gc-cs [21] -4,5695+2.0371i 
-0.09+0.00011 
Gc-pPr opsoed -5.0382 Gc_pso [31] -3.4175+2.1215i 
Gc-ga [20] -4.9100 Gc-3rd pr opsoed -23.8219+0.0000i; -3.6415+14.4301i 
Gc-pso [20] -4.8800 Gc-sa [21], [31] -26.7093+0.0000i; -3.5354+13.9157i 
Gc-csa [21] -4.3670 Gc-gr[21], [31] -26.6997+0.0000i; -3.5402+13.9077i 
Gc-pso1 [30] -4.6300 Gc-mr[21], [31] -19.5655+0.0000i; -4.2123+14.6057i 
Gc-Propsoed1 -4.9182+1.4597i Gc-cs[21], [31] -25.1041+0.0000i; -3.6830+14.0186i 


7. CONCLUSION 

The bycycorobot model balancing is achieved using the reduced order controller methodology. The 
sixth order controller is pre-specified in the literature. The controller when comes closed-loop with the 
system did not perform well, while in the complete closed-loop system, performance gets improved, while 
the order of the system in complete closed-loop increases to sixteen orders. Therefore, a reduced order 
controller is projected and successfully obtained using the HHO in the first, second, and third order. The 
result in a complete closed-loop is compared with the result of the higher-order controller, reduced order 
controller. The proposed controller not only balances the bycycorobot but also gives better performance in 
the first order, second order, and equivalent performance in third order as its performance is compared with 
the reduced order available in the literature. 
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